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ABSTRACT 

An  investigation  is  made  of  the  drift  motion  of  small  particles 
under  the  influence  of  acoustic  oscillations.     The  investigation  is  made 
to  determine  if  the  motion  has  a  magnitude  great  enough  to  produce 
significant  changes  on  the  fuel  distribution  in  the  chamber  of  a 
liquid  propellant  rocket  motor.     The  calculations  are  made  for  the 
motion  in  both  a  rectangular  and  a  cylindrical  chamber. 

In  the  rectangular  chamber  the  gas  oscillation  is  restricted 
to  the  fundamental  transverse  mode  and  motion  in  only  one  dimension 
is  considered.     The  particle  drift  velocity,  that  is  the  non -oscillating, 
non-damped  term  in  the  expression  for  the  particle  velocity  is  found 
in  the  solution  of  the  second  order  equation. 

For  particle  motion  in  the  cylindrical  enclosure,  only  gas 
oscillations  in  the  first  transverse  or  sloshing  mode  is  considered 
and  motion  is  restricted  to  a  transverse  plane  of  the  cylinder.     The 
particle  drift  velocity,  again  a  second  order  term,  is  determined. 

The  magnitude  of  the  drift  velocity  is  calculated  using  condi- 
tions found  in  a  liquid  propellant  rocket  combustion  chamber.     Dis- 
tances a  typical  fuel  droplet  would  move  during  its  average  life  time 
are  calculated.     The  distances  are  small  compared  to  the  size  of 
most  rocket  combustion  chambers. 
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I.     INTRODUCTION 

It  is  a  well  known  fact  that  liquid  propellant  rocket  combus- 
tion chambers  often  experience  a  phenomenon  called  unstable  com- 
bustion.    The  gas  properties,   such  as  pressure,  density  and 
temperature,  fluctuate  during  periods  of  unstable  combustion. 
The  fluctuations  are  not  random  and  definite  correlation  exists 
between  the  fluctuations  at  different  positions  in  space  or  different 
instants  in  time.     Characteristically,  the  fluctuations  have  a  fixed 
frequency  and  the  amplitude  of  the  variations  can  grow  to  values 
comparable  to  the  steady  state  value  of  the  property. 

The  regularity  and  the  large  amplitude  of  the  fluctuations 
distinguish  the  property  variations  observed  during  unstable 
combustion  from  the  random  variation  in  gas  properties  observed 
during  "steady"  operation  of  a  combustion  chamber. 

The  observed  frequencies  of  oscillation  vary  from  10  to 
as  high  as  10,000  cycles  per  second  and  depend  on  the  driving 
mechanism  or  phenomenon  and  the  chamber  size.     The  frequen- 
cies can  be  divided  into  two  general  categories:    (1)    a  low  fre- 
quency or  chugging  instability  with  frequencies  in  the  neighborhood 
of  100  cycles  per  second  and  (2)  a  high  frequency  or  screeching 
instability  with  frequencies  in  the  range  of  1000  to  10,000  cycles 
per  second. 

The  low  frequency  phenomena  are  often  caused  by  an 
interaction  between  the  rates  of  propellant  injection,  combustion 
and  rate  of  ejection  of  material  from  the  chamber. 
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In  the  high  frequency  case  the  variations  of  gas  properties, 
during  unstable  combustion,  are  accompanied  by  motion  of  the  gas 
in  the  chamber.  In  general,  the  motion  of  the  gas  and  variations 
of  properties  are  similar  to  the  conditions  which  would  be  gener- 
ated by  a  strong  acoustic  oscillation.     This  observation  suggests 
that  the  large  amplitude  variations  observed  during  high  frequency 
unstable  combustion  are  a  result  of  a  resonance  which  occurs 
when  a  driving  mechanism,  associated  with  the  combustion  phe- 
nomenon, is  in  phase  with  the  natural  acoustic  modes  of  the  com- 
bustion chamber.     For  example,  in  cylindrical  chambers,  the 
oscillations  occur  in  the  radial,  tangential  and  longitudinal  modes 
plus  their  various  combinations.     The  first  transverse  or  sloshing 
mode  is  the  mode  of  oscillation  most  frequently  found  experimen- 
tally    'in  combustion  chambers  with  high  frequency  instability. 

The  present  investigation  will  be  concerned  with  the  high 
frequency  instability.    A  complete  and  satisfactory  answer  to  the 
question  as  to  why  this  instability  occurs  has,  as  yet,  not  been 
found.     In  the  attempt  to  answer  the  above  question,  many  mech- 
anisms have  been  proposed  that  are  based  on  the  interaction  of  the 
pressure  and  temperature  oscillations  with  the  combustion  process. 
The  interaction  depends  on  the  existence  of  a  time  delay  which  is 
essentially  the  delay  time  in  the  response  of  the  burning  rate 
to  a  change  in  temperature  or  pressure. 

Since  the  delay  times  in  combustion  systems  are  very 
sensitive  functions  of  the  fuel-oxidizer  mixture  ratios,  a  systematic 
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variation  of  the  mixture  ratio  in  the  chamber  may  set  up  conditions 
that  are  favorable  for  unstable  burning.     The  suggestion  has  been 
made  that  variations  in  the  fuel -oxidizer  ratio  can  be  produced  by 
the  relative  mean  motion  of  droplets  of  fuel  and  oxidizer  under 
the  influence  of  acoustic  oscillations  in  the  combustion  chamber 
gas. 

This  investigation  is  concerned  with  the  initial  require- 
ments for  the  production  of  the  mean  particle  motion  mentioned 
above.  The  requirements  are  first,  that  a  particle  drift  motion 
does  occur  in  a  chamber  with  acoustic  oscillations  and,  second, 
that  the  magnitude  be  large  enough  to  result  in  the  appreciable 
particle  motion  which  would  be  necessary  to  produce  significant 
changes  in  fuel  composition  in  the  chamber. 

The  following  two  idealized  problems  will  be  solved. 
First,  as  a  simple  example,  particle  motion  in  a  rectangular 
chamber  with  acoustic  oscillations  in  the  fundamental  transverse 
mode  is  considered.     In  the  second,  more  complex  problem  the 
calculation  is  carried  out  assuming  a  cylindrical  enclosure  with 
acoustic  oscillations.     The  small  particle  drift  velocity  induced 
by  the  first  transverse  mode  of  oscillation  is  determined. 

To  calculate  the  particle  drift  velocity  the  equation  for 
the  velocity  of  the  gas  undergoing  acoustic  oscillations  in  the 
chamber  is  combined  with  Stokes  drag  law  to  calculate  the  forces 
acting  on  small  particles.     The  resulting  equation  is  linearized 
by  introducing  a  perturbation  expansion  in  powers  of  a  term 
proportional  to  the  amplitude  of  the  drag  force  acting  on  the 
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particle  as  a  result  of  the  gas  motion.     It  is  necessary  to  go  to 
second  order  terms  in  the  solution  for  the  particle  position  in 
order  to  obtain  a  mean  drift  velocity. 

Representative  values  of  the  drift  velocities  are  calculated 
using  conditions  found  in  a  liquid  propellant  rocket  combustion 
chamber.     Directions  and  distances  the  droplets  would  move,  in 
their  average  life  time,  are  determined.    Finally,  an  upper  limit 
is  obtained  for  a  maximum  variation  in  particle  density  as  a  re- 
sult of  the  drift  velocity. 
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II.     RECTANGULAR  CHAMBER 
Development  of  the  Equations.    Before  proceeding  to  the  more 
complex  and  more  realistic  problem  of  the  investigation  of  par- 
ticle motion  in  a  cylindrical  enclosure  with  acoustic  oscillations 
it  will  be  illuminating  to  look  at  particle  motion  in  a  rectangular 
enclosure  containing  acoustic  waves.     The  rectangular  case  will 
illustrate  the  important  features  of  the  problem  in  a  straight- 
forward manner  since  motion  in  just  one  dimension  need  be  con- 
sidered. 

The  method  of  solution  for  the  particle  drift  velocity  will 
be  as  outlined  in  the  introduction. 

In  solving  the  problem  for  the  case  of  the  rectangular 
chamber  the  following  idealized  assumptions  will  be  made: 

1.  The  fluid  medium  is  considered  to  be  at  rest  and  in  a  uni- 
form state  except  for  weak  isoenergetic  fluctuations  in  velocity 
and  state  properties. 

2.  The  small  particles  are  considered  to  be  spherical  and  to 
experience  Stokes  law  forces. 

3.  The  particles  have  negligible  effect  upon  the  gas  motion  in 
the  chamber. 

Consider  a  rectangular  chamber,  as  shown  below,  with 
acoustic  oscillations  in  the  fundamental  transverse  mode. 


°*l ']_ 


%-■£ 
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Gas  motion  is  restricted  to  oscillations  in  the  ru,  direction.     The 
linearised  partial  differential  equation  for  the  gas  pressure  dis- 
turbance in  the  chamber,  under  the  conditions  considered  here, 
is  a  simple  wave  equation  of  the  form: 

The  well-known  solution  to  equation  2. 1,  restricting  the  solution 
to  real  values,  is: 


p   -    P  sin  ^fr  cos  cot 

lor. 


(2.2) 


C 

Here,    /o    is  the  amplitude  of  pressure  fluctuations  of  the  gas  in 
the  chamber,    C    is  the  local  velocity  of  sound,   60    is  27?    times 
the  frequency,  V  ,  and  "t  is  the  time.    CO     must  be  ^jL'ffi     to  sat- 
isfy the  boundary  condition  that  the  gas  velocity  and,  hence,  the 
pressure  gradient  vanish  at  the  chamber  walls  CV  ~-X  ).    For  the 
first  transverse  mode  considered  here  rn    is  1.     The  velocity  and 
pressure  fluctuations  of  the  gas  are  related  by  the  following  equa- 
tion: 

p1^      -     -VP     -     -^  (2.3) 

The  velocity  of  the  gas  in  the  chamber  is  given  by  integration 
of  2. 3  as: 

U    =    --fir  CoS  ^  Sin  cot      .  (2.4) 
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Equation  2.  4  shows  that  the  gas  velocity  is  zero  at  the  chamber 
walls  and  a  maximum  at  the  center  of  the  chamber.     The  ampli- 
tude of  the  gas  velocity  oscillations  is  proportional  to  cos  JOfe   . 

If  the  location  of  the  particle  is  given  by  f   ,  and  if  asaump* 
tion  2  holds,  then  the  drag  force  on  the  particle  is  given  by  Stokes 
law  and  is: 

(>ir/jcr   (U~-j[)        <  (2.5) 

The  equation  of  motion  for  the  particle  is  obtained  by  equating  the 
inertial  forces  and  drag  force,  and  is: 

styl^I      =    I*/**-  (U-  -#)  (2.6) 

It  is  important  to  note  that  the  gas  velocity  must  be  evaluated  at 
the  position  of  the  particle,  that  is  U  =  U(f.t)    •     If  assumption 
3  holds,  the  gas  speed,  U  ,  can  be  eliminated  by  substituting  equa- 
tion 2.4  into  2.6;  this  gives: 

mn£L    -    -   klAL^-R  cos—i-  Sin  cot  -  &ifa(rAL  (2.7) 

This  result  can  be  simplified  by  the  use  of  the  following  dimen- 
sionless  variables: 

Z   =JT      and  ^T=    ,||    =f±      _  (2.8, 

Introducing  the  dimensionless  variables  in  equation  2.7  gives: 

dJ£    -t  Is'cLZ   =   e'cos  %Z   sin  airy     .  (2.9) 
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Herc,  the  constants  are  defined  as: 


Kb  ?±m*j£     s^    e'«  - ^ ^7 ^ 


(2.10) 


The  solution  to  equation  2.  9  will  give  the  position  of  the  small 
particles,  tE.  ,  at  any  time,  U  ,  while  they  are  under  the  influ- 
ence of  the  oscillating  gas  in  the  chamber. 

Note  that  equation  2.  9  is  non-linear  due  to  the  appearance 
of  2?    in  the  trigonometric  function  on  the  right  hand  side  of  the 
equation.     Because  of  this  nonlinearity,  an  approximate  method 
of  solution  of  the  equation  must  be  used.     The  approximation  used 
depends  on  the  fact  that  the  term   £     can  be  held  to  small  values 
by  limiting  the  magnitude  of  the  viscosity,  or  other  flow  proper- 
ties  appearing  in   £    .     Note  that  this  restriction  on  the  calculation 
is  independent  of  that  required  to  obtain  a  linear  acoustic  equa- 
tion (2. 1)  in  the  first  place.     Thus,  when  the  expansion  in  terms 
of   £    is  continued  to  terms  of  higher  order  than  the  first,  there 
is  no  obligation  to  extend  the  approximation  in  the  acoustic  equa- 
tion to  higher  powers  of  pressure  or  velocity  amplitude.     By  re- 
stricting  the  analysis  to  arbitrarily  small  values  of    £    ,  it  is 
possible  to  linearize  equation  2.9  by  introducing  a  perturbation 
expansion  of  the  form: 

Substituting  the  expansion  into  equation  2.  9  and  separating  the 
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equation  into  equations  for  the  various  orders  of    £    gives: 

iff  +  ^'jt  =  °  (212> 


<£*±    -hk'djL   -  cos%?0  sin  2^  (2.13) 


<&  *-  /<^    -  -  g  2,  sin  %  2o  Sin  slit  T    .  (2. 14) 

In  making  the  above  substitution,  cos  Jk^  is  expanded  as 

For  a  particle  at  rest  at  3=  0   and  at  the  point   2     »  the 
appropriate  initial  conditions  are: 

2e(o)  =   I  j     jp-(o)  =   0 

Z,  (o)    =    &  rt>)     -  o  \  <2' 13) 

Za^)     -  jp  (o)    =    0    . 

Solutions  to  equations  2. 12,  2. 13,  and  2.  14  contain  terms 
which  are  constant,  oscillatory,  exponential  and  linear  functions 
of  the  time.     The  purely  oscillatory  solutions  to  the  above  equa- 
tions are  not  of  interest  since  no  net  transport  of  particles  is 
involved  and  the  solutions  involving  decaying  exponentials  have 
a  negligible  effect  on  the  drift  velocity  after  a  few  oscillations. 
Hence,  the  only  parts  of  the  solution  that  are  of  interest  and  that 
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result  in  a  particle  drift  velocity  are  those  solutions  which  are 
proportional  to  time . 

Looking  at  equations  2.12  and  2,13,  it  is  seen  that  the  so- 
lutions to  these  equations  are  of  the  form  of  constants,  decaying 
exponentials,  and  sinusoidal  terms.     However,  equation  2.14 
has  Z,  •  multiplying  the  right  hand  side  so  that  the  product  of  the 
trigonometric  solutions  of    Z,  with  {s/VjaTrTJ  will  produce  a 
constant  term.     The  solution  for    ZA  will  then  contain  the  first 
and  presumably  most  important  term  proportional  to  time. 

In  the  following  paragraphs  the  solutions  for   zf0  ,    Z, 
and     j^_    will  be  found.     The  chief  aim  of  the  analysis  is  to  ob- 
tain the  single  term  which  is  proportional  to  time  so  that  the 
drift  velocity,  as  defined  here,  can  be  evaluated. 
Solutions  to  the  Equations.     By  inspection  it  is  obvious  that  the 
solution  to  equation  2. 12,  with  the  given  boundary  conditions, 
is  2e  =    Z  =  constant. 

Substituting  Z0-Z     in  equation  2.  13  leads  to: 

i  li     -h  k'd^L    ~    Cos  £  Z  sin  airy  (2.16) 

JO*  &  ** 

The  solution  to  the  homogeneous  equation  is; 

z,  =  C,  +  CAe~       t  (2.17) 

The  trial  particular  solution  is: 

z?  =  C,  Sin  airy  +  C^  cos  any  (2.18) 
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Substituting  equation  2.  13  into  2.  16  gives: 

-  ^7rsCz  siriaird"  -   4  Traces  air 3"  +  K^/rCz  cos  sir^J 


(2.19) 


-k'ar  C    Sin   3.7r^f     - 


COS  ■%  £    Sfn  ^ir^y 


Equating  the  coefficients  of  like  terms  gives  the  following  values 
for  the  constants: 


c,  - 


CoS%  Z 
-4ra-  *'* 


(2.20) 


The  solution  for     ?,    can  now  be  written  as: 

Z,  -|C,  +  C,  C    *K   °°SJ*   5/wOV  4^z  COSttJ 


4ir  -K 


■kTw-Zir' 


(2.21) 


The  initial  conditions  that  2,  (<?)  =    — —7(0)  -  0      are  used  to  eval- 
uate  the  constants   C,     and    Qo    •     ^he  following  expression  is 
obtained  for  the  complete  solution  to  equation  2.13: 


H,= 


nr 


*    £ 


/<  C*S  ^ 


K'^TT+S'Tr3  -4?rV-  K3 


%2 


+   cos*±  Sjnairy  +    Ji-^s  ^    cos  air? 


(2.22) 


The  differential  equation  for    ^    can  be  found  by  sub- 
stituting 2.  22  and  i-0--d.      into  equation  2. 14: 


i  +  *'#  =[-*'***(££?£  %feHj-=^ 


%  Sin%2[-^r-rJ&r)e  sin  air 


5        H^V-*'3 


r 


-f^aF^®S)^/«r3f 


(2.23) 
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The  above  equation  may  be  simplified  slightly  by  using  the  follow- 
ing trigonometric  identities: 

Sin  %E  cos  $2      -    -^SinirZ 

Sm  Zrr^  cosJttT  =    ^sin^rrZT 
(sinairT)S=    ^  -  ->  cos  4?r  CT 
Equation  2.  23  reduces  to: 


,2 


&  _ 


_  / 

(F  V(5'J  Sin.  Q-ir^y   -  F  C~      sin,^7fj 


-J- Sr  ~  %  cos+Try  +  E'smi-TrT    I  (2.25) 


Here, 

'        irs'in  TTZ  '         k'smTt-Z 

F  '=    7r^s''n  r?L-  r  ' -      k'sin  ^^_ 

The  solution  to  the  homogeneous  equation  is  again: 

3,  =  C,  ,Ca  e-  K'*. 

The  following  trial  solution  will  be  used  to  find  the  particular  solu- 
tion to  equation  2.  25. 

Z2   ~  (Cj^y  +  C^  sin  STrCT  -/-  Cs  cos  sury 

_  t'zr  _fr  <2-26> 

-+~Ck£  sin  sir  If +■    C?  e        cos  Stt^T 


■+-  Cy    Oos^iry    -f  Ccf    Sin   ^-if^f  I 


-13- 
Substituting  2.  26  into  2.  25,  equating  coefficients  of  like  terms, 
and  solving  for  the  constants  gives: 


,^/ 


(2.27) 


Mr*-3 v     S/n4-rJ  -hi ~T$~7r~3    sinz-rrJ 

Here  it  is  unnecessary  to  evaluate  the  constants  C,   and    C^  , 
since     C,    is  a  constant  and  will  disappear  when  the  derivative 
with  respect  to  time  is  taken  to  find  the  particle  velocity  and    Ca_ 
is  multiplying  a  decaying  exponential  which  will  not  contribute 
to  the  steady  state  particle  drift  velocity.     The  non-oscillating, 
non-decaying  term,  a  second  order  term,  is    Z^  ~/ — ^T~3l~~^\  J  CT 
Using  equation  2.11,  the  non -oscillating  and  non-decaying  part 
of  the  expression  for  the  particle  position  can  be  written  as: 

z  =  z  +  e*  *?'*?&  ?  +  -  -  -  <2-28> 

Or,  using  2.3,  equation  2.  28  can  be  written  in  terms  of  the  origi- 
nal variables  as: 


?K 


'(**>+«y 


.14- 
Now  taking  the  derivative  of  2.  29  with  respect  to  time,  the  mean 
particle  drift  velocity  is  obtained: 


!  =  -  -  * 


S'n~jP 


Remember  here  that  £' -       ~  ^aRJ  „d|<=3*2^. 

Equation  2.  30  gives  the  steady  state  drift  velocity  of  the 
particle  which  is  at  a  point    f    in  the  chamber.     The  drift  term 
is  second  order  in  £     and  depends  sinusoidally  on  position.     The 
greatest  drift  velocity  occurs  at  an  initial  particle  position  of 

The  fact  that  the  drift  velocity  is  second  order  can  be  jus- 
tified by  the  following  argument.     The  drift  velocity  exists  because 
of  the  gradient  in  the  amplitude  of  the  gas  oscillations.     Consider 
the  forces  acting  on  a  particle  due  to  the  difference  between  par- 
ticle and  gas  velocity.     The  gas  is  oscillating  at  a  fixed  frequency 
but  with  an  amplitude  that  varies  spatially.     Thus,  during  a 
single  oscillation,  the  particle  moves  from  a  region  in  which  the 
amplitude  of  gas  motion  is  large  to  a  region  in  which  the  amplitude 
of  gas  motion  is  smaller.     During    that  part  of  an  oscillation  when 
the  particle  is  near  its  maximum  excursion  in  the  direction  of  in- 
creasing amplitude  of  gas  motion  forces  act  to  accelerate  the 
particle  toward  the  direction  of  decreasing  amplitude.     Similarly, 
when  the  particle  is  on  the  opposite  side  of  its  mean  position, 
forces  act  to  accelerate  the  particle  in  the  direction  of  increasing 
amplitude  of  gas  motion.     However,  since  the  force  acting  on  the 
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particle  is  proportional  to  the  difference  between  the  gas  speed 
and  the  particle  speed,  the  forces  acting  on  the  particle  are  higher 
in  the  regions  of  higher  amplitude  of  gas  motion.     Hence,  over  a 
number  of  oscillations  of  the  particle,  a  net  force  acts  to  accel- 
erate the  particles  toward  the  direction  of  decreasing  amplitude 
of  gas  oscillations.     The  net  driving  force  depends  on  the  differ- 
ence between  the  forces  acting  on  the  particle  on  either  side  of 
the  mean  particle  position  and  thus  depends  on  the  mean  difference 
in  the  gas  speed  on  either  side  of  the  particle  mean  position. 
Hence,  the  driving  force  and  the  resulting  drift  velocity  is  pro- 
portional to  the  product  of  the  gradient  in  amplitude  of  gas  motion 
and  the  amplitude  of  particle  motion.     Since  both  of  these  quan- 
tities are  first  order,  the  resulting  particle  drift  speed  is  second 
order. 

The  gradient  in  the  amplitude  of  gas  oscillation  is  pro- 
portional to  sin  —-2    and  the  amplitude  of  the  particle  oscillation 
is  proportional  to  cos-3-rz?     ,  e.g.   see  equation  2.22.     Thus,  the 
maximum  drift  speed  would  be  expected  at  the  value  of    Z     which 
maximizes  (cos  JrZ)(S'n?'  zf)         or  ?  =   ±Je£  . 
Discussion  of  Solutions.     It  is  of  interest  to  determine  the  man- 
ner in  which  this  drift  depends  upon  the  particle  size.     To  do 
this,  it  is  convenient  to  use  the  fact  that    £    can  be  written  as  a 
function  of     /<»(£  =  ° £    )   •     Therefore,     £    can  be  elimi- 

nated from  equation  2.  30  to  give  the  following  expression  for  the 
drift  velocity: 


/?a     K 
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X 


To  make  the  dependence  upon  particle  radius  explicit,  write  the 

particle  mass  as   rm.=  Q  %1T <s-       ;  then  J\  -  ^zr~z — 5.  •     I»  these 
terms  the  drift  velocity  is  given  by: 


As  the  particle  radius,  c  -  $  gets  large  it  is  evident  from  equation 
2.  32  that  the  drift  velocity  approaches  zero.     If  equation  2.  32  is 
rewritten  as: 

ra     SinJZl  <2.33) 

it  is  evident  that  as  the  particle  radius  gets  small,  the  drift  ve- 
locity again  approaches  zero.     Therefore,  a  maximum  drift 
velocity  exists  for  some  finite  value  of  particle  radius.     This  op- 
timum  radius,    o      ,  resulting  in  maximum  drift  velocity,  can  be 
found  by  examining  the  first  and  second  derivatives  with  respect  to 
c   .     The  result  of  the  calculation  gives: 


*    =  W_  (2.34) 

Substituting  the  optimum  particle  radius,  2.  34,  into  the 
drift  velocity  expression,  2.  32,  gives  the  following  expression 
for  the  maximum  drift  velocity: 
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* 


SL,  )  sin  iX 


In  terms  of  the  chamber  pressure,    fc    »  *ae  maximum  drift  is: 

Note  that  while  the  optimum  particle  radius  is  a  function  of 
the  chamber  dimension,  x  »  *ne  maximum  drift  velocity,  rj",  is  in- 
dependent of  the  chamber  size.     Therefore,  the  drift  motion  is  less 
important  as  the  chamber  size  increases. 

Note  also  that  the  maximum  drift  velocity  is  independent  of 
droplet  properties  and  only  depends  on  the  gas  properties  through 
the  speed  of  sound   C   ,  specific  heat  ratio   y    and  pressure  ratio 

t  " 

The  ratio  of  drift  velocity  for  arbitrary  particle  radius  to 
the  drift  velocity  for  optimum  particle  radius  is  simply  obtained 
from  equations  2.  32  and  2.  34  and  is  given  by 


or 


(2.37) 


The  drift  speed  drops  rapidly  for  particle  radius  different  from 
optimum  radius.     For  example,  for  6~"  -  J:  e-  or2r    the  ratio 
^£*  is  about  y  . 
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m.     CYLINDRICAL  CHAMBER 
Development  of  the  Equations .     The  techniques  used  in  Section  II 
are  now  applied  to  the  solution  of  the  problem  of  particle  motion 
in  a  cylindrical  geometry.     The  problems  are  quite  similar  except 
that  with  the  cylindrical  geometry,  (A)  motion  in  two  directions 
must  be  considered  and  (B)  centrifugal  forces  may  be  important. 
To  investigate  the  motion  of  small  particles  in  a  cylindri- 
cal enclosure  under  the  influence  of  acoustic  oscillations,  the 
same  assumptions  that  were  made  for  the  rectangular  enclosure 
will  be  used;  namely: 

1.  The  fluid  medium  is  considered  to  be  at  rest  and  in  a  uniform 
state  except  for  weak  isoenergetic  fluctuations  in  velocity  and 
state  properties. 

2.  The  small  particles  are  spherical  and  experience  Stokes  law 
forces. 

3.  The  particles  have  negligible  effect  upon  the  gas  motion  in 
the  chamber. 

For  the  conditions  given  in  assumption  1,  it  can  be  shown 
(3)  that  the  gas  motion  is  governed  by  the  wave  equation  which  is 
written  below  in  cylindrical  coordinates  for  the  pressure  distur- 
bance: 

Since  we  are  interested  in  finding  a  particle  drift  velocity 
in  a  cross  sectional  plane  of  the  cylinder,  in  analogy  to  the  simple 
transverse  waves  considered  in  the  rectangular  chamber,  only 
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oscillations  independent  of    2    will  be  considered  in  this  analysis. 
The  solution  to  equation  3. 1  for  a  circular  cylinder  with  standing 
acoustic  oscillations  is  given  by: 

P=    f?  Tni™?)  X  (ft?)  e"  .  (3.2, 

OX 
Here,  as  in  the  previous  problem,  the  frequency,  ))    ,  is  jf^      , 

and,  in  this  case  the  angular  velocity,  (x\  ,  is     J^J^-~-     .     The 
c^>rn  n  arise  from  the  condition  of  zero  radial  velocity  at  the  cyl- 
inder walls.     The  subscripts  an  and  m  specify  the  modes  of  vibra- 
tion in  the  (f    and  /I  directions.     At  the  walls  cKm™   *s  a  solution 
to 

___        __  0         ^  (3  3) 

For  each  value  oi/m   there  exists  an  infinite  number  of  solutions 
to  equation  3.  3  since  3  ,„    is  an  oscillatory  function  of  its  argument. 
Characteristic  values  of  c%v^  are  listed  in  reference  3. 

For  oscillations  restricted  to  a  transverse  plane  of  the 
cylinder  the  simplest  cases  are  the  first  radial  mode  \pn-=  c  ,  ^  =  /  , 
0^^=  1.  22),  and  the  first  transverse  or  sloshing  mode  {rm-itoi-o  , 
cX^  .586).     In  investigating  the  possible  particle  drift  velocities 
induced  by  the  acoustic  oscillations  only  the  effect  of  the  sloshing 
mode  will  be  considered  since  it  is  this  mode  which  is  most  often 
observed  in  liquid  and  solid  propellant  rocket  systems  experiencing 
high  frequency  instability. 

With  the  restrictions  to  oscillations  independent  of   2     and 
to  real  solutions,  equation  3.  2  can  be  written: 
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p  =    p  cos  (f  3",  (1^)  cos  srrvT 


(3.4) 


The  coordinate  system  used  in  this  problem  is  shown  in  the  follow- 
ing diagram: 


In  this  case: 


VP  =    tM  (sinq>)  3,  (p^)  cos  Sr»t  j 


+  L 


\Bicos<f)(co,S^t)(^)[  X(23£)-  ZJ  (Mrf)]  (3.5, 


An  expression  for  the  gas  velocity  can  be  obtained  from  equation 
3.  4  in  a  manner  analogous  to  that  used  in  the  case  of  the  rectangu- 
lar chamber.     The  result  of  the  calculation  is: 

R 


U  =  tv  [=>£$.  s!n  f  J'<  (^— )  sin-  S7rvT] 


R 


^[-^k  *»<Ps/n***tfe(2±±)_  ^(?r)}l         (3.6) 


Fig.   1  shows  the  magnitude  and  direction  of  the  gas  velocity, 


as  a  function  of 


a 


H, 


,  for  various  points  in  the  cross  section  of 


the  cylinder.     The  calculations   are    made  for  a  time  such  that 
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sin  9.W"0t  -  1.    Since  the  cross  section  is   symmetrical,  only  the 

first  two  quadrants  are  shown  in  the  plot.     It  can  he  seen  in  Fig. 

o  o 

1  that  there  is  a  gas  velocity  node  at    (f  equals  0    and  180    while 

the  pressure  node  occurs  at   LP  *  90    and  270   . 

Using  assumption  2  the  drag  force  on  the  particles  is  given 

by: 


6  7T/<  6-  (U  -nrj  (3.7) 

Here,  as  before,     (J  is  the  gas  speed  and   nr  is  the  particle  speed. 
The  particle  position  is  given  by    i^  J\   ,  the  particle  velocity  by 
oVt  —  -f    dik     •+•  ~l(n  a,  dJL  »  an<^  ^e  Particle  acceleration  by 

f£   -   t   fj%    -»/d<t'fl    ,   -T  [3J/ld£     ,   A  dXW   7 

7t   ~  "LJt*       (it J  J  +  V/*5F57  +  *-&}  ■ 

The  equation  of  motion  for  the  particle  can  be  obtained  by  equating 
the  particle  inertial  forces  to  the  drag  force  obtained  from  Stokes 
law.     The  result  is: 


t         Id 


%-*ffi] +*,"»[*&%**&] 


(3.8) 


+t  {^t cos  * sin  *™Tfe  w- T*  (^  -#}]. 

In  the  above  equation  the  same  symbols  are  used  for  gas  position 
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as  for  particle  position  since  the  gas  velocity  must  be  evaluated 
at  the  position  of  the  particle  (by  assumption  3).     Equation  3.8 
can  be  separated  into  two  more  useful  equations  corresponding  to 
the    i^    and     6^    directions: 
I  ro    component 

6\  ZiryUe-    a/a    _   a  /Jcfif'     - 

At*      +       nm  ~Jt  l  \dt  J        " 

3f?ri£VcL  T  LJ"\    a    /    Jil    a    /J 

L  to  component 


(3.9) 


*£$  +*&  +  *%?**  1^  S,'nl?Sin  ^^  T,(^ i     (3<  10) 


As  in  the  case  of  the  rectangular  chamber  it  will  be  conven- 
ient to    nondimensionalize  equations  3.9  and  3.10;  the  following 
dimensionless  variables  are  used: 

Introducing  the  dimensionless  variables  into  equation  3.  9  and  3.10 
gives: 
1*.    component 

j^£    +  K  dr  t  ^rJ  (3.12) 


m  £,3  — 

k 

cf    component 

-  (3.13) 

£  SM  0  5/n  ^^  T  ("*«,.*)  J   ^ 

Here,  £  =  3yts<r°R  ^3   ^_  i_2^*I£      are  parameters  of  the  prob- 
lem  and,  except  for  constants,  are  the  same  as  the    &    and    K   in 
the  rectangular  chamber  problem. 

Equations  3. 12  and  3.13  are  non -linear  due  to  the  appear- 
ance of    z?   in  the  argument  of  the  Bessel  functions  and  are  analo- 
gous to  equation  2.9  in  the  case  of  the  rectangular  chamber.    Since 
£    is  proportional  to    PQ  t/U ,  and  other  flow  properties,  an  argu- 
ment similar  to  that  made  about  equation  2.  9  is  valid.     The  follow- 
ing perturbation  expansion  in  powers  of    o    is  introduced  to  linearise 
equations  3. 12  and  3. 13. 

z(y)  =  ?0  (?)  +  e  z,  (r)  +  ea4  (?)  +  --  -  P.  *4> 

<pM    =    <Po  (?)    +  £  <fr  (?)  +  f<t±  (?)  +  -  —  (3. 15) 

In  substituting  the  above  expansions  in  equations  3. 12  and  3. 13, 

sin  p   is  expanded  as    S'ln  <f  =  (sin  fe)(l-  -|  [$,  +  £& T  + ) 

■+■  (cos$)(£  eft   +£  <f>a )  and  cos  <p     is  expanded  as 

cos^(cos^)(i-^x+e^f^.--)  -  (smtyJltQ+eX---)  . 

The  Bessel  functions  are  expanded  in  a  Taylor  series  about  the  point 
2C  to  get: 
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To  (irot,9i)    =     To  (  t<*„*0  -  e  ?,  Tf*l0  T,  (r«Kl0  H„)  -£a£,  7r<*,e  ^C^Z,) 


-  e3?* 


(r^f  td  (7r«,020)  ♦  ^§^  j;r^,.2.)] 


3;  (-»-=<,„  2)  = 


Z  (*"<.  ?J   +"  £  Z,  "<*<•  ^  (^Z.)  -  ?p  Jj  (r«„^) 


+  eaZa  ^,0  t,  (tw„a)  -  ap>  ^  (^2^  ^S%-fcf  3,(ir<x,.Z 

The  results  of  substituting  3. 14  and  3. 15  into  equations 
3. 12  and  3. 13  and  separating  the  results  into  orders  of  £     are: 
L  /l    component 


(3.16) 


dO7"         *d?"  °  dVdT 


oCT 


^   i 
/J 


(3.17) 
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6 


age  sin  zjn^y  r3Z  ^  c^5  ^  -j;  ^  Ztf) 


c 


_  22, 


COS  fc  ^a  (Tr~„Z.)  +  <p,  Sin  <P0  f  yo  ( v-k  2J  -  ^(Y^)]    (3. 18) 


/  (0   component 


aZo5r5r  +  ^^  +***$  =  °  (3,19) 


^JyJy  ^tjyjr     Uzf>dy  sy  +?°  jy*  +2ttjr* 


■*-<rt, 


+  Ze/C^   +  Z 


Ki!  ^   =  [jA  VJ)  s/nfc  sz/i  ^fa^  ]  O.20) 


dJdj        jyjy         dy  jy         dTd^r        jy  jy      ^dy^f 

*  <(>,  cos  p0 1,  CTr^ftij  S//i  ^aCr      . 

If  a  particle  is  considered  to  be  at  rest  at  >-  -  C    and  at  a 

point  ■-  ,    -      ,  then  the  appropriate  initial  conditions  are: 

2,  Co)  =  ^*(o)    =    (f  (o)   =  £k(o)   =,   c  <3-22> 

jy  '  dy 
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It  is  noted  again  that  purely  constant,  damped  or  oscillatory 
solutions  will  not  be  of  interest  since  no  net  transport  of  particles 
is  involved. 

The  form  of  equations  3. 16  to  3.  21  can  be  compared  with 
equations  2. 12  to  2.14.    Since    (b (    ,  and     ?,     multiply  the  right  side 
of  the  second  order  equations  it  is  expected  that  any  terms  in  the 
solution  which  are  proportional  to  time  will  again  appear  in  the  sec- 
ond order  equations  3. 13  and  3.  21. 

The  expectation  that  the  drift  velocity  will  be  second  order 
in    £  ,  specifically  for  the  radial  motion  of  the  particles,  also  fol- 
lows from  the  physical  arguments  given  on  page  14.     However, 
for  the  geometry  considered  here,  centrifugal  forces  also  act  on 
the  particle.     Since  the  centrifugal  forces  are  proportional  to  the 
square  of  the  particle  velocity  and  the  particle  velocity  is  first 
order  in   £   ,  the  drift  which  results  from  centrifugal  forces  is 
also  a  second  order  term. 

Solutions  to  the  Equations.     A  solution  to  the  zero  order  equations 
3. 16  and  3. 19  is    z?0  =  constant  =    ?    and    <P0  =  constant  =  (J)    .     The 
constant    £    is  closely  related  to  a  drag  force  and  can  be  called 
an  interaction  parameter  so  that  to  zero  order  in   £    there  is  no 
drag  or  interaction  between  the  particles  and  the  gas  in  the  cham- 
ber with  the  result  that  the  particle  position  is  fixed.     That  is  to 
zero  order  in  £    ,  the  particles  are  not  moved  by  the  gas,  but  are 
at  a  fixed  position.     Using  these  simple  zero  order  solutions,  the 
solutions  for  ?,  ,    Z^   ,    $,    ,  and  0^  can  be  found  by  the  same 
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method  used  in  solving  equations  2. 13  and  2. 14  in  the  rectangular 
chamber  problem.     The  complete  solutions  to  equations  3. 17,   3. 13, 
3.  20  and  3.  21  are  discussed  in  more  detail  in  Appendix  A. 

The  first  and  presumably  most  important  non-oscillatory 
and  non-damped  terms  appear  in  the  solutions  to  the  second  order 
equations  3. 13  and  3.  21.     These  non -oscillatory  and  non-decaying 
solutions  are: 


2*  = 


&  = 


\0    -  ^(FE  +  DG)  I  y 


A 


(3.23) 


Here, 


L 

E 


-    :x  (r^,0g)  sin  $ 


=    Tf£*  cos}  jra  (Tr<*,0  *)  - 10  (iro^z  )J 

H      =       Sin  $  [r«*,o  T0  (ir^f )  -^r  31  (ir«*teZ)] 

j      -       cos  $  HT  0r«^0!) 

Differentiation  of  equation  3.  23  with  respect  to  time  will 
give  the  particle  drift  velocity.     The  objection  might  be  made  that 
the  decaying  exponential  terms  in  the  solutions  to  the  first  order 
equations  may  contribute  appreciably  to  the  particle  motion.    A 
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c ale ulation  to  determine  the  magnitude  of  the  first  order  decaying 
exponential  terms,  compared  with  the  steady  state  drift  velocity, 
is  shown  in  Appendix  B.     In  general,  these  terms  appear  to  be  neg- 
ligible for  reasonable  values  of  the  parameters. 

Moreover,  the  terms  which  decay  exponentially  in  time 
arise  from  certain  assumed  initial  conditions  of  the  particle  posi- 
tion in  a  certain  phase  of  the  motion.    Since  there  is  no  physical 
way  in  which  one  initial  condition  or  phase  may  be  preferred  over 
another,  these  can  occur  at  random  and  hence  the  decaying  terms 
can  not  contribute  to  the  mean  motion  of  large  masses  of  particles. 
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IV.     RESULTS  AND  DISCUSSION 

The  Drift  Velocity.     Ignoring  the  oscillatory  terms  and  decaying 

exponentials  the  drift  velocity  of  the  particles  can  now  be  written 

down  using  the  solutions  listed  in  3.  23.     Remember  that  the  velocity 

of  the  particles  is      /V    --    t  17    -f  i/i^ 

at  T      c/t 

or  in  the  dimensionless  expanded  form  is: 

(4.1) 
The    (-i    component  of  the  drift  velocity  can  be  written  as: 

in  terms  of  the  original  variables.     Note  that  [_-££>  -  (Ft  +D(r>)] 
is  a  constant  for  any  given  /[    ,  m    and  cylinder  radius  and  does 
not  depend  on  the  particle  dimensions.     L-et  [^D  -  (F£~+DG)j  =  $ 
then  the  /I    component  of  drift  velocity  is  more  conveniently 
written  as 

"* =  zliKi^fi     .  (4-3> 

To  see  how  the  Si,    component  of  particle  velocity  varies  with  par- 
ticle size  the  expressions  for   £    and     ^    can  be  substituted  into 
nf^     1  this  substitution  gives: 
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^   - 


(4.4) 


4y<n'^3 


Or,  since  the  mass  of  the  particle  equal.3    P  JZ  If  <r      ,  4.  3  can 
be  written  as 

As    6"-^oo   in  equation  4.5,  /7£— >  0    and  as    &-*— >•  o  »/^->C! 
therefore,  as  in  the  case  of  the  rectangular  chamber,  there  is  an 
optimum  particle  radius,    cr     ,  which  results  in  a  maximum  <o£   . 
This  optimum  radius  is  again  found  by  examining  the  first  and 
second  derivatives  of  the  velocity  with  respect  to  6" .     The  result 
of  the  calculation  gives: 


The  same  arguments  can  also  be  made  for  the    £p  component. 
Thus,  the    ($>   component  of  particle  velocity  is 

in  terms  of  the  original  variables.    Again,  ^a.  [EH+Dl)  is  constant 
for  a  given  A    ,    0     and  cylinder  radius.     If  ^q.  (B H + Dl J  5    <£  ; 
then  the    (P   component  of  drift  velocity  can  be  written  as 

~£  C&,    .i-rrr/  <48> 


*      \3LK[lc*  +  pg*Tj  J 
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Since  the  expression  for  the  /I  component  of  drift  velocity  (equa- 
tion 4.  3)  differs  from  the  if  component  of  drift  velocity  (equation 
4.8)  only  by  a  multiplying  constant,  it  is  evident  that  the  optimum 
radius  for  maximum   (V^  will  again  be  <r ■   -  J-^f&yZ.  • 

Since  the  optimum  particle  radius  is  the  same  for  the  A> 
and    (f  components  of  drift  velocity,  a  particle  of  radius    (T      will 
indeed  experience  a  maximum  resultant  drift  velocity.     Note  that, 
since   V^=  — F    i»  *ne  case  of  the  rectangular  chamber,  the  optimum 
radius  given  in  4.  6  agrees  exactly  with  the  optimum  radius  found 
in  equation  2.  34  for  particle  motion  in  a  rectangular  chamber. 

Again,  in  analogy  with  the  rectangular  chamber  problem, 
it  will  be  interesting  to  see  whether  or  not  the  maximum  particle 
drift  velocity  is  independent  of  the  cylinder  radius,  a  .     The  par- 
ticle drift  velocity  can  be  written  as: 


Iiere   2>     and    $   contain  the  cylinder  radius  in  the  ratio  /q.    «"*<* 

for  a  given    -1z    they  are  independent  of  the  radius.     It  remains  to 

be  shown  that  £  ,  K    and  (tq^    *r©  independent  of  the  cylinder 

radius.     Since  \)  =  c^'°-      ,  it  is  obvious  that   V    varies  as  the  in- 

^  di- 
verse of  the  radius,  and  that       %JOQ\     does  not  vary  with  the 

radius.    Since  the  optimum  particle  radius  varies  as   CLa     and  the 

mass  of  an  optimum  sized  particle  varies  as  CL     ,  it  is  evident  that 

£  =    3Xc^Jg.     and  |<=   £%^g~g    are  both  independent  of  the  cylinder 

radius.     Therefore,  while  the  optimum  particle  radius  is  a  function 

of  the  cylinder  size,  the  magnitude  of  the  maximum  particle  drift 
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velocity  is  not.     This  result  is  in  agreement  with  the  result  ob- 
tained for  the  rectangular  chamber  problem. 

Another  comparison  with  the  rectangular  chamber  problem 
can  now  be  made.     If  the  optimum  particle  radius,  equation  4.6,  is 
substituted  into  the  particle  drift  velocity  expression,  equation  4.9, 
the  following  is  obtained: 

**=  Tfa^P  (%)*[ 1*  ~  f»  Sl  ■         H-10> 

For  cp  =  0  *  the  constant    ^     is    Q    and  the  particle  motion  is  en- 
tirely in  the  radial  direction.     Fig.    1  shows  that  the  gaa  velocity 

o 

in  the  cylinder  is  also  entirely  in  the  radial  direction  when   (f    is 

Since  motion  was  restricted  to  one  dimension  in  the  rectangular  geometry 

problem,  a  comparison  of  the  maximum  particle  drift  velocities  is 

o 

most  meaningful  if  compared  for  a    (D  of  0   in  the  cylindrical  problem. 

0 

When    (f  -  O  is  substituted  into  equation  4.  10,  the  equation  reduces  to: 

«- *  C  ff  (if fe'7^  -^^fa%IZt»&-Zt*$  <4-  »*> 

For  comparison,  equation  2.36  is  rewritten  here: 


r-A^b]  SinZl  (2.36) 


x 


The  expressions  are  the  same  except  for  the  terms  on  the  right  re- 
sulting from  the  geometry  difference  in  the  two  problems. 

It  is  also  noted  that  the  ratio  of  particle  drift  velocity  for 
arbitrary  radius  to  the  drift  velocity  for  particles  of  optimum  radius 
varies  in  exactly  the  same  manner  as  given  in  equation  2.  37  for  the 
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rectangular  chamber  problem. 

Magnitude  of  the  Drift  Velocity  and  Discussion.     Thus  far,  with  the 
assumptions  listed  in  Section  II  and  the  restrictions  mentioned  for 
£,  the  expressions  for  particle  drift  velocity  are  perfectly  general. 
The  expression  for  particle  drift  velocity,  equation  4.9,  will  now 
be  used  to  calculate  the  drift  velocity  of  small  liquid  droplets  in  a 
liquid  propellant  rocket  combustion  chamber.    Although  the  fluid 
medium  in  the  combustion  chamber  is  not  at  rest,  the  iviach  number 
of  the  gas  relative  to  the  particle  is  usually  low  near  the  injector. 
Stokes  law  drag  will  hold  if  the  droplet  shape  remains  nearly  spher- 
ical and  if  the  Reynolds  number  of  the  flow  is  less  than  one  which, 
for  reasonable  values  of  the  parameters,  means  that  the  particle 
radius  must  be  less  than  ten  microns.     The  results  of  these  drift 
velocity  calculations  will  not  be  exact,  but  will  give  an  indication 
of  the  mean  motion  the  droplets  experience  in  a  combustion  cham- 
ber. 

The  following  properties,  for  the  droplets  and  gas  in  the 
chamber  will  be  assumed: 
For  the  particle: 

£p     =  2. 1  slug/ft.  3 

The  particle  radius  will  be  picked  to 

maximize  the  drift  velocity  and  for 

the  values  listed  here: 

_*"  .5 

6       =  2.  75  x  10       ft.   =  8.4  microns 

For  the  gas  in  the  chamber: 

C  =  2360  ft. /sec. 

/a  =  1.54x  10       slug/sec.  ft. 

P     s  4.  lh  x    10"4  aluo/ft.  3 
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T  =  2500°  K 

\>  =  692  1/scc. 
For  the  chamber: 

a  =  1  ft. 
In  substituting  the  above  values  in  the  equations  the  form  of  £ 
was  changed  to  £  =  1%~fr-H    by  substituting  V  ^    for     C 
Here,     f£     and   P  are  the  mean  chamber  pressure  and  density. 

In  calculating  the  magnitudes  of  the  drift  velocities,  values  of 
A   of  .1,   .5,  and  .  9  ft.  were  chosen  for  the  initial  particle  posi- 
tion, and  since  the  drift  velocities  are  completely  symmetrical, 
values  of  (f    in  the  first  two  quadrants  of  the  cylinder  cross  sec- 
tion were  considered.     In  particular,  angles  of  0    ,   22.5    ,  45    , 
67.  5    ,90    ,  and  the  corresponding  angles  in  the  second  quadrant 
were  used.     Table  I  shows  the  results  of  the  drift  velocity  calcu- 
lations.    In  both  Table  I  and  Fig.   2  the  magnitude  of  the  particle 
drift  velocity  is  divided  by     \/P)        •     Since  the  amplitude  of 
pressure  fluctuations  in  an  unstable  rocket  chamber  is  a  parameter 
which  varies  quite  drastically,  it  is  convenient  to  multiply  the  values 
listed  in  Table  I  and  Fig.   2  by  the  square  of  any  desired  ratio  of 
pressure  amplitude  to  mean  chamber  pressure  to  determine  the 
actual  drift  velocity. 

Fig.   2  shows  the  directions  and  magnitudes  of  the  drift 
velocities  in  the  first  two  quadrants  of  the  cylinder  cross  section. 
It  is  seen  that  the  drift  near  the  wall  of  the  cylinder  is  greatest  at 
Cp  =  90    which  is  to  be  expected  since  at  this  value  of  (fi    there 
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is  a  pressure  node  and,  as  can  be  seen  in  Fig.    1,  the  gas  velocity 
is  entirely  in  the    (p   direction.     The  drift  velocities  are  generally 
in  the  radial  direction,  due  to  the  centrifugal  force  exerted  on  the 
particles  by  the  oscillating  gas.     The  drift  velocities  are  greatest 
when  the  initial  particle  position  is  at  the  mid-radius  of  the  cylinder. 
Since  the  drift  velocity  is  a  second  order  phenomenon,  it  is  difficult 
to  find  any  correlation  between  the  drift  direction  and  magnitude, 
and  the  pressure  and  velocity  distribution  of  the  oscillating  gas  in 
the  chamber  except  for  that  mentioned  above. 

To  obtain  a  range  of  values  for  the  magnitude  of  the  drift 

p 
velocity,  a  value  of  .  316  is  chosen  for      /G     \  that  is,  a  pressure 

amplitude  amounting  to  31.  6  per  cent  of  the  mean  chamber  pres- 
sure.    Then    (  Vp  )       is  0. 1  and  with  the  values  shown  in  Table 
I,  the  particle  drift  velocity  varies  from  .  39  to  4  feet  per  second. 

As  mentioned  earlier,  the  magnitude  of  the  first  order  de- 
caying exponential  terms  in  the  solution  of  the  cylindrical  chamber 
problem  are  investigated  in  Appendix  B.     For  this  numerical 
example,  the  results  show  that  after  just  two  oscillations  the  mag- 
nitude of  the  velocity  contribution  of  the  decaying  exponential 
terms  is  less  than  one -one  hundredth  of  the  magnitude  of  the  second 
order  drift  velocity. 

It  is  also  interesting  to  determine  how  far  the  particles 
move  in  the  combustion  chamber.     In  .01  seconds,  which  is  close 
to  an  upper  limit  for  liquid  droplet  life  times  in  a  rocket  chamber 
(4),  the  droplets  can  move  from  .  05  to  .43  inches  depending  on  their 
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initial  position  in  the  cylinder.     This  is  a  very  small  distance  when 
compared  with  the  3ize  of  most  rocket  chambers. 

In  comparison,  the  drift  velocity  in  the  rectangular  chamber 
will  reach  a  maximum  of  17.  4  feet  per  second  at  an  initial  position 
in  the  center  of  the  chamber.     In  .  01  seconds  the  particles  could 
move  2  inches  which  is  considerably  greater  than  in  the  cylindrical 
chamber,  but  still  not  significant  when  compared  with  chamber  size 
or  wave  length  of  the  oscillation. 

Variation  in  Particle  Density.     Although  it  is  not  the  aim  of  this 
present  work  to  calculate  in  detail  the  variation  in  fuel -oxidizer 
ratio  which  results  from  the  relative  drifts  of  fuel  and  oxidizer  par- 
ticles, it  is  of  interest  to  calculate  a  few  representative  values  for 
the  change  of  particle  density  with  time.     The  variation  in  particle 
density, m  ,  can  be  computed  directly  from  the  continuity  equation 
for  nonsteady  flow6.     The  equation,  in  vector  form,  is 


^m 


+  S7-(mnr)   ~  0     .  (4.12) 


J* 

nr  is  the  particle  drift  velocity  and  m  =m (t, fif  Cf)  is  the  droplet  density 
in  number  per  unit  volume.     If  the  problem  is  restricted  to  the  cal- 
culation of  variations  in  density  over  short  periods,  then  the  varia- 
tion of  /7l  with  time  may  be  approximated  by 

Also,  for  short  periods,  the  effects  of  evaporation  of  the  particles 
may  be  neglected  and  hence  the  particle  radius  v/ill  be  constant. 
Vvith  the  above  restrictions  equation  4. 12  reduces  to 
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A/?%?     ^    -(v^)^"f     •  (4.13) 

Values  of      sfae    are  calculated  for  the  chamber  conditions  given 
in  the  previous  section.     The  calculation  is  carried  out  for  the  cyl- 
indrical  chamber  at    six.   =0.1  and   if  =  0  and  for  the  rectangular 
chamber  at   )/fc    *  0. 1.    A  time  of  about  1/5  of  the  particle  life  time 
of  0.  01  seconds  is  assumed  and  a  radius  equal  to  the  optimum  value 
for  particle  motion  is  used  in  the  calculation. 

A  3  per  cent  reduction  in  particle  concentration  is  obtained 
for  the  cylindrical  chamber  and  a  10  per  cent  reduction  for  the 
rectangular  chamber.     These  values  of  density  variation  are  the 
maximum  which  can  be  obtained  for  the  chamber  conditions  used. 
For  the  cylindrical  chamber,  the  changes  in  particle  density  appear 
to  be  too  small  to  cause  a  large  change  in  oxidizer -fuel  ratio  and 
in  fact  the  variations  of  the  order  of  5  per  cent  are  much  less  than 
the  variations  which  occur  because  of  the  injector  design. 
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V.     CONCLUSIONS 

It  has  been  demonstrated  that  small  particles  in  an  oscil- 
lating gas  field  experience  a  steady  state  drift  velocity  in  both 
rectangular  and  cylindrical  enclosures.     This  result  is  based  on 
the  assumptions  that  (1)  the  fluid  medium  is  at  rest  and  in  a  uni- 
form state  except  for  weak  isoenergetic  fluctuations  in  velocity  and 
state  properties,  (2)  the  particles  are  spherical  and  experience 
Stokes  law  forces,  and  (3)  the  particles  have  negligible  effect  upon 
the  fluid  motion. 

The  particle  drift  velocity  is  a  maximum  for  certain  opti- 
mum sized  particles.     The  optimum  particle  radius  is  given  by 
0-  _     /Ji£_     where    P    is  the  frequency  of  acoustic  oscillations 
in  the  combustion  chamber. 

The  maximum  particle  drift  velocity  is  independent  of  the 
chamber  size  although  the  optimum  particle  radius  is  a  function 
of  the  acoustic  frequency  and  hence  of  the  chamber  dimension. 

The  magnitude  of  the  particle  drift  velocity  in  a  typical 
liquid  propellant  rocket  chamber  varies  from  .  4  to  4  feet  per  sec- 
ond and  in  an  average  particle  or  droplet  life  time  the  actual 
movement  will  amount  to  just  a  few  tenths  of  an  inch. 

It  is  concluded  that,  for  rocket  chambers  of  reasonable 
size,  the  drift  movement  of  particles  has  a  small,  if  not  negli- 
gible, effect  on  the  distribution  of  fuel  or  oxidizer  in  rocket  com- 
bustion chambers. 
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TABLS  I 
PARTICLE  DRIFT  VELOCITIES  IN  THE  CROSS  SECTION  OF  A 
CYLINDRICAL  ENCLOSURE 


7L 

'a 

(degrees) 

(ft/ sec) 

(ft/aec) 

(ft/aec) 

0 

11.8 

0 

11.8 

22.5 

10.6 

-.3 

10.6 

45 

7.5 

-.4 

7.6 

67.5 

4.5 

-.2 

4.5 

90 

3.9 

0 

3.9 

.5 

0 

37.4 

0 

37.4 

.5 

22.5 

34.4 

-5.3 

34.8 

.5 

45 

27.1 

-3.7 

27.4 

.5 

67.5 

19.7 

-2.6 

19.9 

.5 

90 

16.7 

0 

16.7 

.9 

0 

10.9 

0 

10.9 

.9 

22.5 

12.2 

-6.9 

14.0 

.9 

45 

15.4 

-9.8 

18.2 

.9 

67.5 

18.5 

-6.9 

19.8 

.9 

90 

19.9 

0 

19.9 
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APPENDIX  A 
SOLUTIONS  TO  EQUATIONS  3.  16  TO  3„  21 
In  this  appendix  solutions  are  obtained  for  equations  3.  16 
to  3,21  which  are  developed  in  Chapter  III, 

If  a  particle  is  considered  at  rest  at  3^=  o  at  a  point   Z  ,    (p  , 
then  the  initial  conditions  for  the  equations  are: 


z.w«  2  ,  <p0(o)  =  <|>    i|s(^  =  jAp(o)  = 


JT 


jr 


0 


(i) 


Zero  Order  Equations.  The  zero  order  equations,  3.  16  and  3.  19 
are: 


0 


(2) 


(3) 


z?0  =  constant  =  2    and    cf*©  =  constant  =    <J>   satisfy  equations  2  and 
3  and  will  be  used  as  the  solution. 

First  Order  Equations.    The  first  order  equations  reduce  to  the  iol- 
lowing  upon  substituting    Z    and     (J)    into  3.  17  and  3.20: 

t  y    component: 


2ai^  +  f  ^  J|,  ,Pr  (To,g)  sm  ?  s//j  ^p^- 


(4) 
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L  x  component: 


3- 


^  s/n^sJ'cos  <f  [j;C7rcv?)  -Titrate)] 


(5) 


The  solution  for  <pt    is  obtained  first  and  then  the  solution  for    i? ( 
is  found. 


If 


Ds 


s/'/i  <f> 


(6) 


the    Im     component  equation  may  be  rewritten  as 


c/*<fc      -4-    l<f  ill      = 


D  s/ai 


c 


(7) 


The  solution  to  7  is  determined  by  the  usual  technique  of  finding  the 
complementary  and  a  particular  solution  and  adding  the  two,     The 
constants  appearing  in  these  solutions  are  then  eliminated  by  use 
of  the  initial  conditions. 


The  complementary  solution  to  equation  7  is: 

,-kO'] 


TIC 


c,  +  c 


For  the  particular  solution  try: 


<5  7TPa  J 


(8) 


(9) 


In  this  appendix,   the  subscript   C     will  pertain  to  the  comple- 
mentary solution  of  the  equation  and  the  subscript    p    will  pertain  to 
the  particular  solution. 


T,p 


4>  = 


The  first  and  second  derivatives  of  equation  9  are; 


C3  Co,  il 


G4     S'" 


(10) 


-  (^f  C3  Sin  V^  -  (^/C4  cos  & 


(11) 


L 
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Substituting  10  and  11  into  equation  7  results  in: 

C  3  C  C  4-  C 

Equating  coefficients  of  like  terms  gives: 

-(£***/ C3     -K^Br   Q      =     D, 
and 


.   (iz) 


(13) 
(14) 


The  solution  of  equations  13  and  14  is: 


C    --  ^ 


c  = 


-  ;<  D 


k*  +  (ig2f-     ,        ^  -  affi2rf^(5|Si/j      . 


The  complete  solution  for    (J)      is  then: 


(15) 


4>  = 


C,    +  C^C  -r 


-D  sin.  —  ^ — 


/<£  cos 


WiT 


C 


[k<+  (»p.y] 


(16) 


To  determine    C,     and     C     the  initial  conditions  are  used, 


<P,   (o)    =»  O   = 


4>%;  a  0  = 


C,   +  C^    -  - 


KD 


¥*[**+&Ff] 


-KG, 


D 


c. 


k*+  ( 


a£va\s 


(17) 


(18) 


The  solution  of  17  and  18  is 
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c.  = 


D 


ixyi 


KD 


k[kC(ZF)'\  **%*[&  m*] 


The  complete  solution  for    d>     is  then: 


(19) 


«  = 


D^ 


-f  - 


K  D 


3TV0L. 


_    kX 


k[?^f]   ^Ovmr]    w+i'^y] 


D 

5  in 

y 

ka  + 

err 

_     KD  cos 


£V£<4. 

c 


fK.+  (^Ha-/J 


(20) 


A  similar  process  to  that  carried  out  above  leads  to  the  solution  of 
r,     from  equation  5* 


If        E  = 


\"JA  (lro<?)cos  $  -  T6  (7r*2-jGqS  <£  ][  (21) 


then  equation  5  reduces  to: 


E  sin  aTrvQ-7 

u  c 


(22) 


Since  this  is  exactly  the  form  of  equation  7  and  since  the  initial  condi- 
tions on    <p     and  j?(   are  identical,   the  solution  can  be  written  down 
immediately  by  analogy  and  it  is  given  by: 


-47- 


r~ 


2,     = 


£T   2-rry>a- 


K  E 


t  -3—  e 


E  sin  ^^ 


/<  E    COS    -3 — 


<3W4 


[k%(^r]_ 


Second  Order  Equations       The  second  order    iK    component  equation, 
3.  18,   reduces  to: 


Z,  [ffii) \  T,  (Tor?)  COS  4> 


L  <z 
1 


when  zero  order  solutions  are  substituted. 


and 


F-= 


6=-fe   T„  (  tt«J)  s/«  ?  -  T&  ^  <  ***)  s"7  * 


then  equation  24  may  be  expressed  as 


24) 


(25) 
(26) 


J  3" 


I  +  K  4$&  -  Z(j$f=  \zt  F  5/«  ^p-r  f-  4>  fe  s/«  a»v«r 


d4>> 


(27) 


The  values  for      ~JJ£l,    ,     7      ,   and  d>      can  be  obtained  from  the  solu- 
tions  given  in  equations  20  and  23.     Thus,   from  equation  20: 


-48 


r 


J4\ 


D 


(Jjrp:)£-kr-   ^Z^D  cos**?*?   +  kOSmii^J 


^  -frV^N-*- 


k^  (^) 


-  2^2JP  Sin  *Jg*? cos  ¥g±?   t  k*DzSin*2X&2 


~1 


[k^    (i-^)2] 


21*- 


(28) 


(29) 


Substituting  29,    20  and  23  into  equation  27  gives: 


r?Da(^)2e-ikT-    2  2D*(te^Je-kZc^Jr^ 


2  JTt?^\2'l^ 


[k*+  m&fj 


o 


■f 


•f    * 


[ fc ^  ( dgS/j  k  L  K *  +  ( iZp)*] 


k[k^(^-^y] 


■49- 


k 


£  F  OOb  ~Z        i"i  ~3 4_    K  D(v  Sin  *~^r- 


*[K  *+(*%**)*] 


0  6  Sin***?*? 


K06  o 


0,  €lp±J  SjK  affKUy 


[**+(*&*)*] 


(30) 


Equation  30  can  be  simplified  by  using  the  following  trigonometric 
identities: 

Sin  — £        cos      g  _    /2  bm      £ 


S/**-*?^ 

c^ 


s  (4  --i«5-t^3-) 


(31) 


C0£ 


a  ;rP*^      = 


The  results  of  substituting  the  above  identities  into  equation  30  is: 


*  -2k* 


a  .  Kt 


<^    +  ^   -  |gDV4^)V        -  ^D^^fe     ^-^ 


c/r: 


^r 


L 


[K*+     (i^rfT 


lk*+(U^?] 


+  (FB  +D6)(^)sia*2p?    ^  (Ft+DG)kSin&l^Lr 


k[k>+  (2™*)*] 


Lv+(*%*yj 


-50- 


(F£  +D6)  k/^    Sin 


c 


37rV<x 

c 


[K't  C-^Y] 


(32) 


Equation  32  is  the  linear  second  order  differential  equation  which 
must  be  solvedo     The  solution  is  obtained  in  the  same  manner  as 
was  used  in  the  previous  section  to  obtain  values  of  the  first  order 
functions    <p     and    £?,   . 
The  solution  of  the  complementary  function  is: 


and  the  trial  particular  solution  is: 


(33) 


-t  o7  y°s  —£  +  O&  s'^  — 3; —     +  C<?  cos  — s —    +  L>l0sm   — 3— 


(34) 


The  first  and  second  derivatives  of    i?       are: 


?'    - 

•^ 


~KCLe       sin  ± 


-t 


■    rc    -oes  r— ,         _   j^ 


(35) 


2*  = 


-51- 


kT. 


-KT        t7Tv*"^ 


+  K*JO!«  C;Z     s/«  ^Jg^  -  m**fcs  C       cos  * 


_rr 


tr^1«^-a-w  ,<Cle-K3c,sW- -K^ctcvS-%? 


-(&&)%&-     sihZg^  -(igs-fCjC****^  -{ilpYc^i 


1*     c  .'^  4-r^ti  J 


(L> 


-{*M±fC,  os  argj  -  (tUgifC,  u* 


27r>><<'J 


(36) 


Substituting  35  and  36  into  equation  32,   and  equating  coefficients  of 
like  term  results  in  the  following  evaluation  of  constants: 


kc,- 


2  k.  C4  - 


k^K-^,3 


,  \  <a  flVa 


(37) 


(38) 


(39) 


(40) 


(42) 


\     c,     J    u\ 


27ryz. 


(43) 


-52. 


2  vy±  i<  r 


"(-f^ 


c 


10 


(ht:  +D&)lA  -j-  ((•  £+DG)  k 


(44) 


Solving  equations  37  to  44  for  the  constants  gives: 


C. 


a 


J   I  2TrV4_\  a- 


1 


'(Fi+UG)(*VLl 


(45) 
(46) 


c7  = 


^(^V^+^V^       -*Jp(F£+D*) 


\[*9+  (^frPTfk'+itg*)^     [k\  (^jji^pi^^i 


k5 


t-jiiEil^l 


) 


c 


([*•%  (^)'f^ [^r^py]  [k'ti^TTU'+e^F?] 


(47) 


(48) 


cf  = 


-  (Ff  ^P6) 


n 


(ztts)YK"+(^?")7 


C,„- 


-    (F£  +D  G) 


IK  (^P)l^  i2-^)'] 


(49) 


The  complete  solution  for    z?j  can  now  be  written  as: 


r 


2,  = 


C,    +  Cj  £ 


+  l**cf-(Fe+DGJ]c:**c<K  ZTS  +  ff^M)gJe-xs/,^ 


/>*+  (2?^ 


k[k3+  (3jyj±.y] 


-53- 


(50) 


To  determine    C,      and    C     the  initial  conditions  are  substituted  into  50: 


\M  -  0   -- 


C,  f  C,    + 


ZQUB&zy*   ^    a%D*_  (fz+qg) 


jk 2im jf ( m&yf  lk*+i ^y] * 


-/-      ^    (  F£+Pg) 


(Ft  #■  oe; 


] 


r?*^  f^^  (svwvtf-  ]  3 


i^A^fFTi^c-tu      (&?*•)*  [k'+  m&yj 


(51) 


F,  (o)   -  £ 


_K[jj>D>-(F^D6)]     +    (FB  +  DG)  (*&*)*+   [%*(*?>fjCDl-%k'c*] 

i**+  m^yj*      K[\d+(^nT  kMvpyyiwipf] 


+  ^  K  (f£+D&)  -    I p£  +v&) 


Solving  for    C,     and  Q     >   results  in: 


(52) 


r 


c.  = 


L 


-54- 


vhh  +  DG). 


k'CK»+(^PJaJ 


(53) 


c.= 


[4  Da  -jJrfr-E+oOl  -f  goa(a^g;Vfag[f-(pr»D<p)3 
(fi^;aj  *a[/ca+  (^p)3Ja  f  **+  f^g^j * 


_  kaI>a  + 


+   (FE4D0)       -  _J_D M2^!     -    [??Pa-  (FE-fSG)]^     (FE+P6) 
-  [-%  (a^)V  i?r"a  *-azfcalfapJ    +  ^f '  (FE+P6  ) 


-i 


_     ^  (Fc  +  £6^ 


5rm 

c 


r^rfj^w-rrJ 


(54) 


The  complete  solution  for     I?      is  obtained  by  substitution  of  equations 
53  and  54  into  50  and  is: 


2,= 


,  I  <_ 


-  (FE+D6)(H?S)*.    fez^V-fefc'o']    -%(Ff+P6) 


-55- 


(Ft  +  DZ)  FPa(^r^    [3*P*-   (FE  +  D6)] 


K/i   (f£+C 


_ -t_ (Fg+DG)         -f  L%. P- -£(F£*K>}£ 


-W 


•.  •,  1  i 


-    k^^)2]*         f*M*?W  k*lK'+(*%±yj 


t 


o  \Z     -  3  A.  T 


2. 


+ 


^  K  (F<=  +0J0)  SM  t^  J     _  (F*  +oO  ^- ^  "       _   (F£  fD6)5-> 


[&(^f]^[&(+Tti   erfi^(*Y*rl    k(3T)L^(^T]  J 


(55) 


-56- 

The  second  order    Ly    component  equation  can  be  treated 
in  a  similar  manner.     Thus,    substitution  of  zero  order  solutions  in 
3,  21  leads  to: 


If 


/-/  =    sw  $  L^  Jo  (nw-2)  -£  T,  l^zf ) J 


and 


(57) 


J     =    [j,f7T^zf)oCS^] 


(58) 


then  equation  56  can  be  expressed  as: 


iSigLih  +lai^  +2f^d^  +K?a&+-aKZzM  = 


jr  jr 


o* 


JT; 


o/J 


Xj 


2,  H  s/h  3M*3-    -^  j  s/n  awarj 


(59) 


Terms  involving     z?,    and    (£,     are  evaluated  by  use  of  equations  20  and 
23  and  are  given  by: 


M<  - 


o 


(60) 


(61) 


57- 


kJ 


+  jD(^)KEe-      sm^l^   ±(Z!^fDECof*Zi£ 


[ka+  (^)2j 


Jli 


-  2  (*£*».)  DKE  s\nilO&?&!>zJry>±X    -f-  /<5CE   S/'n3  3  try  a  y 


.27* 


[V*  (2^)^ 


(62) 


X63) 


Lk**m&F] 


2,  A  - 


_kct 


- k ~ED  (UP*-)  eT  "  V  (£&*)  k£C  Sin  ^  S  J  +  "2*  /<  j^   o03  ££>i«  ? 

[k'+CJipy]*' 


[k'+^Y-T] 


a1r^  fl^ff  ^  /  atr^-n-'  ;  * 


[K*  +  (jL^y] 


(nP;  m  s/«  apliozupz   _  fag^/cto  ^a  i»i? 


[ka4(qp)a] 


(64) 


■58. 


2  i$± 


&D(£pjAe-^-  (£&*)* ED  cos  1^+  KB D^^  sin  ^^AJ 


,21   ^ 


+  *SD  (Vp)e-ky  -   (i^)KEDccs^^-i-  k^b  sin  ^a  s 


Kik'-t  m&y] 


•>  1   ■»- 


fcj- . 


l_j^     } -  <•  __C C 


-KtOsin&p?     -KED(^p)e-K:cosZ2p? 


[K*  +  {qp*f] 


*Vi>±\21* 


^UM^'Vj 


27  * 


+  ^KH  D  cos*^^?   -  KXB  C  sin  *™l?c*s  &£& 


*2p    J"ka    +(£^jaj 


Substituting  62,    64  and  65  into  equation  59  gives: 


-3DB  (*?**)  V  2k *++o(Vg±) 2£e 


(65) 


strra-y 


CCS 


lKjA^/yj 


-W. 


-59- 


■f 


AK£D(^)2e-Kr-2E(^)30s>«*-^-*(^?)\k0cn^ 


2  <[*'+&¥*)'] 


ZlDl*-)2!3" 


c 


2  iJP  [kU  C-^Yf 


2  k  Ik'tim^)*]* 


-2keo(^)&sin  3-*P? -H^fED  sk 


*[k' ■>  (*£&)''.] 


mi^y       -  ,.2,  r-  z.:^^,.-^.,  ^';^' 


+*(«?£)  fee  D  sw  ^ws  ^  J    -  3  k3£P<^>     c^^P 


2   IV-M^)7 


*^)/>m^p)T 


pa-y 


-ak'eo  £¥*)&-** '+a(^P)ic*£Dc*s2n>*? 'sfeo  sin  II^J 


2(zgz)[k>+(zg*r]* 


•  60. 


+2\<*ED{^r~)z~K^z^    -t  afceoKVaJe"** s'm  &g& 


21^ 


?*£kM*Ss)'J 


*1* 


*  [k2 +(**?*)*] 


+ak'BD(*3¥s)er  Ka>s  ^  y  -±PjP)**e 0  <*s* &&  J 


+2&D  Sin B2p* cc 5 2™S2    +  {tH+Dl)^  Sin&gSZ 


*WI>M*p*)aJa 


kL**+  ('-¥-?] 


+  (EH+DJ)  k  SI*  ^  r  -  (Bh  +01)^  e    *\)n  *» 


*-&*  [k*+ m&y] 


i<[kj+  (zg*yj 


-p^J)s//ig^  -(EH+DZ)kc*S&g^_  Sjn.  &g      ' 


[*'+(*&*?] 


^T/aiks-t(i^yj 


(66) 


Using  trigonometric  identities,   equation  66  is  simplified  to  the 
following: 


JT2 


-61- 
f  3  (2Z&)  K£D  SintlpJ    -  2  (£*£*)*£ D  sin  *Zg±? 


-M*¥*)KED  Sin  &p&+2  m^)3BD  C'KTSm  &£ 


-2k3EDsin  Zip?  +K3ED  sin  i^  +  (^.  Dj)(2*p)s)n*g*T 
t  (EH+DI)  k  sin  ^y  -  (Eh+D^i^e'^sm  3-^ 


-£  (E  *  +  P3 )  +  £  (EH-hDZ) cc6 *&g*2    _  K^(EH-tDl)  sin  *&& 

tk*7(^)r]  ~3%B[k>+li  "J    (6?) 

The  complementary  solution  is  again: 

K    =[C,     +Ca£.-kJJ      #  (68) 

The  trial  particular  solution  is: 


->C?tos*lp?+  C?  sih*g&  +  Ci  cos  i2p?+  Cl9  sin  Up?  " 


(69) 


The  derivatives  of    &      are: 


ti,sfo  "«i;e"J'-lc(Ve-|t,w  a3ga7-  w«c.e-^^ 


-62- 


-  X  CL  e '  KySin  ^D^2-  +  £_p-  cb  e  '  ^Cos  ifJ.  £2g9>  Q7  Sin  *m3* 


+m&  C9  Cos  &¥?  -  &&*  Cq  Sin  ^^  +-2*ea-  Cio  CoS  22D22J 


-J 


(70) 


ft      r 


€  - 


rip 


4*t4e' 


_i/tj* 


trtve"*^  aga3"+ag»KC»e-*sn  aga? 


+/(*§»  c$  e'^w^-  (affile"  *<Us  ^^2' 


-crfc^     "*:"'''  -  (9^)aC7cos±Q^./;^3cis/ii^p 


\     C       / 


[■   C     > 


{*P'?Ci*s3lg33'  -  (™*fCl0  Sin  V%S3  j 


(71) 


Substituting  the  above  results  into  equation  67  and  equating  coeffi- 
cients of  like  terms  we  get: 


KG: 


(72) 


■3  K  Ua    " 


--DE  f^J 


3-    1 


j?UJ*(*P"t)f_ 


(73) 


tfWr-*3?^* 


-J*EP(*?J)  +£^F/Dt  -(^)(Lb+U) 

z[&[&F)f  zKLkHzspff  k[v%(^)1] 


(74) 


63. 


-(*P)%    "*(*?*)  C*    = 


4-Dlz  pyo-j  x 


(75) 


-^/C7+fc(4?»)C,= 


(76) 


-/<(^rjc7-(^^/c,= 


-  ^(tHHi; 


^[k'H^'Y] 


-t 


Kdt.D 


j  i«? 


3  (*%*)[*>  ti*ysr] 


(77) 


-(^)  C,  +  k(^)c„   = 


o 


(78) 


-fcfWWH^JX  = 


-,   /J  TTV^N3 


-f-  _A£_L 


(f  * ♦  or) df*)     -    afe  *£ o +  H  i e « * ' il 


Kl^H^f]         2(>Jp)[K>+(3{#rT    ^'^(TY]_ 


(79) 


Solving  equations  72  to  79  gives 


c3- 


(80) 


ct-  = 


r-aoE(4?!af  -   (eh+di)   -  att'gp 


zfr'+iupf]1   ik'+(wrr  H^WJ 


(81) 


-64- 


(82) 


C;= 


2  (^)5£D    _     -4  (EH +01)  -^k'ED 


frM^ffO<=H*#J  fk"+p^r/[ic'+«»5»)aj  z^^Kfa'+CWJ 


+  ^ 


%.  (EH+DI) 


-    ICED 


^(^[kw^k*^     ^(^(^[^^Xyf^VWj 


(83) 


c= 


-3k£D(^) 


+  &  (EH^-DI)  -3(^ya)KtD 


^(^[k^^^^^j  (^)f^p¥flfRW¥?  ^(^JTkV^S 


-  k'3ED 


4-  k/i  (fcH+Dl) 


?^%km^-)Tl^t?}    mcvH^fitt'+ppo'] 


(84) 


C,= 


UMeSseSL  t^:;i 


^  i^*i 


-r      a.^££ 


JL  (EH+Pl) -J(  M'tHHl 


(85) 


c,= 


The  complete  solution  for     (T      is  obtained  by  substituting  the  above 
constant  into  equation  69  and  adding  the  complementary  solution. 


(86) 


-65 


<fe  = 


-2DE  maC'kJc«  ^^  -akVoe-  ^c.5^-r-fEH.Ol)e-K^^-r 


^[K^m^r] 


Oa  +  (22p*)a] 


n; 


*1* 


zk  [k\(d%>sy] 


k  [k*+  (*g*-a)*] 


-fC^P    COS   +*£**  -JkEO^fsin+Zp?       +%(EH  +  Dl)s,n*Lp? 

Zlwqpfl'lk^PlFfl  2(¥?)L*x+(^)T{kH^$i  (^[W-TllM'^f] 


(tH-+W/ wi — -r —        —  k   {EH -tOT) cc$-  ^  T^lYytDSn-,; — 


;     3 


[K*+(ag*)a]a        P^/fkM^'j*       hk  [kM*¥*)*j 


•  66- 


(EH+DDksi**^ 


mik'-r^n* 


(87) 


To  determine    C      and    Q     the  initial  conditions  are  used: 


cfyo)  =  o  = 


c,  +  ( 


o 


?k3[K*+(zp)f  *[&+(2S*T]*  lk2+(*g*)2]* 


±a  (*&*)*  ed _->^-  (i  h  +r  a  -_  %_*f£  C 


+    *&  (£H*0*j 


-      A,  ^  D 


+  zed  {zgiy- 


('jp)(^)[w^r]it\e**j]  ^^x^)[KH^)%+m]  ?Om^f)j3 


-r  3KaSD  +  5  **££  -   k*(EH+Dl) 


(88) 


$£<0)  =0  = 


■ku,  -  *  i-cM- 


(£H  +  DI)      +a-DfJ!gz)L-  +  +£  EO 


klK'A^T]     **[**■*{*¥*¥]     *  [K*+(£g±Yj 


oJ 


-t   fe  (kh  +  CI) 


*[pA^x]*[Ka+(±gs?]    [W^Yll&iHpy]    z{k3A3-inTL^(^Y] 


J   3 


-  m**)k*EL  t  %  (*%*)(£ »■+»*)     +  +(*z "  !l*i     — ](89) 

The  constants    G,    and   C/     are  found  from  88  and  89  and  are: 


c*  = 


k*ft*+(ap?)*j      Sfctk^iap*)*]*      *[Ka+(flF*)*J3 


-67 


-oZEDt^r       +i(EH+Dt) 


-3(!lf*)(vr&)£D 


2[k>^jpy]Xe*tp^)x]    [&(>■¥*)'][*'+£&*)*]    2[KaA^)T[K*H*2p?] 


(90) 


c,= 


fitC^f       -   hiEH+CI) 


*os4v^fUMPF$]  [kat(^)3]i^m^rj  f[Kj+r^)!]a[^(i^>j 


-/-  (Ifffa)  K  ^p -    jtb+Cl) 


-±li~£)M_ 


2k*[k:\(2p*y-]     2[k'+(3%*?]3      [icv^ipayj*      2[*V(qp)Ja£|cVff2«J£j 


±   *£(£// *DJ)  ^  ^  /<J£D  -  ^  (EH+DJ) 


f^(if*)lj[/tV(*t*rj    #ftvqp?)3>V(q^J  (^x^W^^flf^^fl 


+  K+ED 


-  2£Q(UpY     -ZK'E  l 


3    13 


^^»jP5»)[^(a$»)a/fjtv^)3      Sfr^PFtfJ3     2  [^(ap?)*  j 


-   jXV/P 


+  k*{£H  +  DI) 


The  solution  for    (f^  is  then  obtained  from  87,    90  and  91  and  is; 


(91) 


-68 


fe- 


k{BHfDl)       -*DE(&Pf     -4-K3£D  +   JED   (4*P±)X 


K3[k2^(i^)1]  ikX^+esp}]*  Zk^^x]3  Zl&P^J'lK^Ctiglfl 


+  3{^)Pg*)ED 


~   ^Uh+H)  +  m&)  k2£0 


Sljc'+i^FJ'fr+pqpr]     [kl+(^%[kLH^'f]   S(^P)IhU^T}\^(^ 


-   (EH+Dl)  ^4-(aJP)3ED      +Ck(^f/af       +•   3 DB  (BJg°) 


:.       i 


[k  ^(^)%KH^f\  *  [f-*(*&*)l3*  *k  Y*  M^ptf    *[*  z+(^p) x] 


+  3(EH+D2)  -iP&*TeO 


+  %>(EH  +  D1) 


[k*+  (i2&rj*    ?lk2+  (^y]iK3+t*g*fl    Oa+  i&pyjik'+PP5?] 


.    +%K'ED  -K^(£HtDJ) 


+  £*£D 


2ftvq^)yJ/^fl^(^^ 


>JED(^)a      -2K'El 


-3k\v 


-*   KJ(h  H  +  11) 


Hf+mpY?  w+^yj3  ^^y\&(2%*ff  {*$*)*[**+ (&&¥? 


~-k(Em-Di)e~K3'     -+■  IDE^^fe-^      +  4-r  el>  e 


j7J 


**[&+(*&*¥]        SfW+ipgsfl*         ?  [**+(*&*)*] 


-pyajk'EDe 


-KJ 


+  (th  fD2)e"k:r        t-fp^/EOe-^ 


?(^)[K^(^ja/fKa+(^rJ   [K^(Ap)aJfKa+^-a/J       ;     KM^M* 


->i  (EH  ♦■01)3' 
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2 K*E  D  e  ~  K'cos  ^^  -  (£  H  -  0 1  )e  *  Jc«  *&*? 


;ia 


-£ (*£*) 2Dt  e   «%in  fg^  +  (EH  +D2  g^ 

-ukB  D(lZp)  c'<JS  in  ^-"^        f-  2  (^-J^j  ^EQ  cos  t  l™^ 


Z  [**<.  {^pY]  °  2LK3+{^?]*h<S  t  (*%**)*] 

-  K  +-F0  cos  ±2g°^  -  2 KED(m^)*5iA  ^Tpy 


+  ^(f//tDl)5,n^  -3(^)^0  Sin  *&*? 


(^Xk'H^fflr+itg*)']        2  [*'+(*&*)*]*  o2*  (*&*)*] 


-K2ED  sin^p^ +-j^  lBH  +  Dt)s>«  *£g±JT 


+  MZgS)kei)s)nUg*3'    -h<2k3£D  i/^rpar  -(BHi-CI )&*■*>)  Sin iM*  J 


-70- 
.    (EHfDX)K  s,^^-Tl 


(92) 


The  terms  proportional  to  time  in  equations  55  and  92  are 
circled  and  are  the  terms  used  to  find  the  particle  drift  velocity. 
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APPENDDC  B 
A  DETERMINATION  OF  THE  MAGNITUDE  OF  THE  FIRST  ORDER 
DECAYING  EXPONENTIAL  TERMS  OF  THE  PARTICLE  VELOCITY 
The  first  order  decaying  exponential  terms  in  the  solution 
for  the  particle  position  in  the  cylindrical  enclosure  are: 


Z.= 


K[KS+  {Sg*tl 


(1) 


Here, 


0     -      HiZ^i)  Sin  $ 
£    =    1C3*  cost  fe(ir*«lJ)   ~^o(rr^/0I)] 

From  equation  4. 1  the  expression  for  the  particle  velocity,  contain- 
ing just  the  first  order  decaying  terms,  can  be  written  as: 


r 


AT    - 


ft  ce^  *?f  ce*  75-1  (2> 


Substituting  the  expressions  for    •£,     and     <£,    into  equation  2  the  fol- 
lowing is  obtained 


nr  = 


For  an  initial  position  of  the  particle  of  i?  -    1  -  *5    ,  =  45    ,  for 

ijfi        =  .  316  and  for  gas  and  particle  properties  evaluated  for  the 


(4) 
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conditions  used  in  the  calculations  for  Table  I,  the  following  ex- 
pression is  obtained  for  the  particle  velocity  due  to  the  damped 
terms: 

nr=k  [7Ue  ;  +  ^(9/.3  e  ) 

For  X  =  .  0023  seconds,  which  is  time  for  approximately  two  os- 
cillations  in  the  chamber,  6  =  Q,        ,  and  the  velocity  amounts 

to: 

nr  =   £  (-•  003)  +  2}  (.oo4)       ^Heo.  (5) 

For  these  same  conditions  the  steady  state  drift  velocity  a3  calcu- 
lated in  Table  I  is: 

&  =  £  (2.1)    +T?  (-.37)       f%ec.  (6) 

Comparing  equations  5  and  6,  it  is  seen  that  after  just  two  oscil- 
lations the  velocity  contribution  of  the  decaying  term  is  negligible 
compared  to  the  drift  velocity. 


